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'_h ' Abstract. We have performed abundance analysis of two slowly rotating, late A-type stars, HD 32115 (HR 1613) 

C/j , and HD 37594 (HR 1940), based on obtained echelle spectra covering the spectral range 4000-9850 A. These 

spectra allowed us to identify an extensive line list for 31 chemical elements, the most complete to date for A-type 
stars. Two approaches to abundance analysis were used, namely a "manual" (interactive) and a semi-automatic 
procedure for comparison of synthetic and observed spectra and equivalent widths. For some elements non-LTE 
^ ■ (NLTE) calculations were carried out and the corresponding corrections have been applied. The abundance pattern 

of HD 32115 was found to be very close to the solar abundance pattern, and thus may be used as an abundance 
standard for chemical composition studies in middle and late A stars. Further, its Ha line profile shows no core- 
to-wing anomaly like that found for cool Ap stars and therefore also may be used as a standard in comparative 
studies of the atmospheric structures of cool, slowly rotating Ap stars. HD 37594 shows a metal deficiency at the 
level of -0.3 dex for most elements and triangle-like cores of spectral lines. This star most probably belongs to the 
5 Set group. 

Key words, stars: atmospheres - stars: abundances - stars: individual: HD 32115, HD 37594 



1. Introduction 

Solar photospheric abundances are usually used as a refer- 
ence for any abundance study. However, it has been shown 
that the mean abundances for some elements obtained for 
large homogeneous groups of stars can differ significantly 
from solar abundances. For example, Luck & Lambert 
(1981, 1985) and later Takeda & Takada-Hidai (1998) ob- 
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tained an oxygen deficiency of 0.3 dex relative to the Sun 
for a large group of A to K supergiants. There are some 
indications that the sun may be slightly iron-rich (+0.09 
dex on average) as compared with solar-like stars of the 
same age (Gustafsson 1998). Extensive abundance stud- 
ies of the superficially normal sharp-lined B- and F-type 
dwarfs were performed by S. Adelman and co-workers (see 
Adelman et al. 2000 and references therein). They found 
that many of the investigated stars show Am-type phe- 
nomena, and only 5 stars with effective temperatures in 
the range 6700-9000 K have abundances close to the solar 
photospheric values. Their study was performed mainly 
in the blue spectral region (up to 5000 A), and therefore 
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abundances of C, N, 0, Na, S, K were absent or rather 
questionable. Varenne & Monier (1999) published abun- 
dances for 48 A- and F-type dwarfs in the Hyades open 
cluster. They have only two stars hotter than 6650 K with 
v sin i < 25 km s _1 . Because of the limited spectral regions 
that were observed abundances for only 11 elements have 
been derived. Half of these abundances were based on 1 
or 2 spectral lines. 

Most sharp-lined stars with temperatures in the range 
6700-8000 K belong to different groups of peculiar stars: 
Ap, roAp, Am, or A Boo groups. They are usually believed 
to have the same or nearly the same atmospheric struc- 
ture as normal stars. However, Cowley et al. (2001) found 
a pronounced core-to-wing anomaly in the Balmcr lines 
of some Ap stars. They could not compare the Ha line 
profiles of these stars with normal solar abundance stars 
because of the lack of the reliable spectroscopic standards 
in this temperature region. We therefore decided to per- 
form an accurate spectroscopic investigation of sharp-lined 
late A- to early F-type stars which are classified as nor- 
mal MS stars, using observations of a wide spectral region 
4000-9850 A. From tabulated rotational velocities (Abt 
& Morrell 1995) we extracted 27 stars classified as normal 
A3-F5-type stars of luminosity classes III-V with v sini< 
25 kms -1 . Among them only 13 stars have metallicities 
in the range \[M}\ <0.15 on the basis of tcmpcraturc- 
gravity-mctallicity calibrations of the observed Stromgren 
colours (a procedure which will be described below). We 
did not consider the remainder of the stars, the majority of 
which have low metallicities [M] <-0.15. In this paper we 
present the results of careful atmospheric parameter and 
abundance determinations for two stars: HD 32115 and 
HD 37594. These stars were classified as A8IV (HD 32115) 
and A8Vs (HD 37594) by Cowley et al. (1969). 

2. Observations and data reduction 

High resolution spectra of the region 4000-9850 A were 
obtained for both stars with the coude-echelle spectrome- 
ter (Musaev et al. 1999) mounted on the 2 m "Zeiss" tele- 
scope at Peak Terskol Observatory near Elbrus mountain 
in Russia. The best resolving power of the spectrograph 
in this operational mode is R=45 000, but because of the 
rather poor seeing we used a wider slit which resulted 
in a reduced spectral resolution of R«30000. A Wright 
Instruments front-side illuminated CCD of 1242 by 1152 
pixels (22 mkm) was used as a detector for simultaneous 
registration of the whole spectral range in a single expo- 
sure. Two spectra were obtained for each star and they 
were coadded after cleaning for cosmic events during the 
reduction procedure. Blucward of 7100 A the consecutive 
orders overlap, while there are gaps of 5 to 55 A between 
the orders for A >7100 A. Typical S/N=150 was estimated 
for the centers of the orders close to the Ha region, while 
S/N decreases to about 100 in the blue and infrared re- 
gions. Note that S/N also decreases from the center to 
both ends of the order. An extra spectrum of HD 32115 
was obtained recently and was used for radial velocity 



measurements, and for equivalent widths measurements 
of a few lines in the infrared spectral region. Table 1 con- 
tains the Heliocentric Julian Dates of the middle of the 
exposure for each pair of observations. The second col- 
umn of Table 1 presents the results of the radial velocity 
measurements (see below) together with the data from the 
literature. 

Spectrum processing was realised with the help of 
a modified version of the PC-based DECH software 
(Galazutdinov 1992). It includes background subtraction, 
echelle vectors extraction from the echelle-images, wave- 
length calibration, continuum rectification, line identifi- 
cation, equivalent widths and radial velocity measure- 
ments. We did not divide the stellar image by a flat-field 
as we found this procedure does not improve the initial 
S/N=100-150 for this particular CCD. 

3. Radial velocities 

Each spectral line chosen for measurement was approxi- 
mated by a Gaussian. We measured 332 and 268 lines in 
the spectra of HD 32115 and HD 37594 respectively. The 
full list of spectral lines in the 4000-9850 A spectral re- 
gion for both stars was extracted from the vald database 
(Kupka et al. 1999). For the final radial velocity we have 
chosen only Fe I lines with the best laboratory wavelengths 
(Nave et al. 1994). Our V T values and their standard devi- 
ations are given in the second column of Table 1 together 
with data collected from the literature. HD 37594 may be 
considered as a single star on the basis of six radial velocity 
measurements, while HD 32115 shows clear V T variations 
and may be a binary system. From the range of the ra- 
dial velocities presented in Table 1 we may conclude that 
our first observations were undertaken near the maximum 
possible separation between the components. Careful syn- 
thesis of the whole spectrum of HD 32115 does not reveal 
any trace of the lines of the secondary star, and therefore 
we conclude that it is most likely relatively faint and thus 
has no significant influence on the primary's spectrum. 

4. Atmospheric parameters 

Stellar parameters derived with different methods are 
given in Table 2. We first derived effective temperatures, 
gravities, metallicities and absolute magnitudes of our 
stars employing Stromgren photometric indices and the 
calibration of Hauck & Mcrmilliod (1998). We used the 
TEMPLOGG procedure (Rogers 1995), which for A3-F0- 
type stars is based on the original calibrations by Moon & 
Dworctsky (1985) with improvements by Napiwotski et al. 
(1993). For the metallicity calibration TEMPLOGG uses a 
relation between the Stromgren metallicity index and the 
stellar metallicity derived by Smalley (1993). The quoted 
formal errors in the stellar parameters from photometry 
account only for the reported accuracy of the photometric 
indices. One can see excellent agreement between abso- 
lute magnitudes obtained here and those obtained from 
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Table 1. Observational log for HD 32115 and HD 37594. Radial velocity measurements based on n lines with standard deviations 
in units of the last digit (in parentheses) are given in the second column. 



HJD-2400000 


K.kms- 1 


n 


Reference 






HD 32115 


52333.2007 


19.15(43) 


67 


this paper 


51929.3035 


42.08(45) 


71 


this paper 


49405.5333 


44.61(43) 




Grenier et al. (1999) 


42675.8313 


43.7 (5 ) 


14 


Nordstrom & Andersen (1985) 


42669.8640 


12.2 (4 ) 


13 


Nordstrom & Andersen (1985) 


42379.7278 


6.8 (6 ) 


14 


Nordstrom & Andersen (1985) 






HD 37594 


51930.2146 


22.52(72) 


69 


this paper 


49736.6142 


22.04(39) 




Grenier et al. (1999) 


42680.8220 


22.4 (3 ) 


15 


Nordstrom & Andersen (1985) 


42675.8928 


24.8 (6 ) 


16 


Nordstrom & Andersen (1985) 


42380.6909 


20.2 (7 ) 


16 


Nordstrom & Andersen (1985) 



Hipparcos parallaxes (ESA 1997). In this case the for- 
mal errors in the surface gravity include both the paral- 
lax uncertainty and an effective temperature uncertainty 
of ±100 K. This provides strong support for the surface 
gravities used in our analysis. Bolometric corrections do 
not exceed 0.06 mag. On evolutionary tracks both stars 
lie close to the ZAMS and have masses slightly more than 
1.5 Mq (Pamyatnykh et al. 1998, Pamyatnykh 2000). The 
age of both stars is estimated at 600-700 Myr. 

To check the influence of the metal abundances on the 
derived atmospheric parameters we calculated a small grid 
of atmospheric models and corresponding colours with the 
abundances typical for our stars (see Sect. 5) and com- 
pared the observed and synthetic colours. Calculations 
were made with the ATLAS9 code (Kurucz 1993) which is 
modified by the inclusion of new line opacities in the form 
of specifically computed ODFs (Piskunov & Kupka 2001) 
and the Canuto-Mazitelli convection treatment (Smalley 
& Kupka 1997). We obtained the best fit to Stromgren 
indices with T off =7250 K and T cff =7170 K for HD 32115 
and HD 37594, respectively. They are given in Table 2 as 
ATLAS9-newODF results. The derived effective tempera- 
tures were also checked by fitting the Ha line profile using 
these models. Special care was taken to achieve a correct 
continuum fit in the cchcllc orders containing the Ha line. 
The typical length of the orders at 6400-6700 A is 90- 
100 A, and we have 2 overlapping orders containing the 
Ha line. The estimated accuracy of our continuum level 
is ±1-2%. The Ha profiles of the models are computed as 
described in Barklem et al. (2000), namely, assuming LTE 
and using the Synth program (Piskunov 1992) with Stark 
broadening calculations from Stehle & Hutcheon (1999) 
and self-broadening from Barklem et al. (2000). Radiative 
broadening and an estimate of the pressure broadening 
by helium are also included. A comparison between the 
observed and calculated Ha line profiles for HD 32115 is 
shown in Fig.l. Plots for HD 37594 show equally satisfac- 
tory fits. No atmospheric model can reasonably fit the core 
of the Ha line in both stars, which is always broader in 



the observations. Atmospheric parameters finally adopted 
for the abundance analysis are given in Table 2. 




6520 6540 6560 6580 6600 

Wavelength [A] 

Fig. 1. A comparison between the observed Ha line profile 
for HD 32115 (thick full line) and those calculated for 3 at- 
mospheric models: 7250_newODF (adopted model - thin full 
line), 7150_newODF (upper dashed line) and 7350_newODF. 



A further important parameter in the abundance anal- 
ysis is microturbulent velocity. It was derived from the el- 
ements with numerous lines and the most accurate atomic 
data. Typically the £ t value as derived from lines of Cr, 
Mn, Fe, Ni, varies from 2.0 to 2.6 kms" 1 in HD 32115, and 
from 2.2 to 2.8 kins" 1 in HD 37594. Ti lines give a sys- 
tematically higher microturbulent velocity in both stars, 
up to 3.2 kms -1 . For abundance calculations we accepted 
Ct=2.3±0.3 kms" 1 for HD 32115, and &=2.5±0.3 kms" 1 
for HD 37594. These values for turbulent velocity fields in 
late A stars agree very well with theoretical predictions 
(Kupka & Montgomery 2002). 
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Table 2. Atmospheric parameters for HD 32115 and HD 37594 with formal error estimates in units of the last digit in parentheses. 



Method 


T e s 


log g 


[M/H] 


M v 




HD 


32115 






Photometry 


7300(40) 


4.21(01) 


0.03(04) 


2.84(02) 


Parallax 




4.24(05) 




2.84(10) 


Ha 


7250(100) 








ATLAs9-newODF 


7250 


4.30 






Adopted 


7250 


4.20 


0.00 






HD 


37594 






Photometry 


7170(50) 


4.21(03) 


-0.15(04) 


2.98(07) 


Parallax 




4.22(05) 




2.90(07) 


Ha 


7100(100) 








ATLAs9-newODF 


7170 


4.25 






Adopted 


7170 


4.20 


-0.30 





5. Abundance analysis 

To begin the abundance analysis synthetic spectra is com- 
puted for the whole observed region with the adopted 
atmospheric parameters and solar abundances. All lines 
that have depths greater than 0.5 % of the continuum 
flux in the synthetic spectrum were extracted from the 
VALD database (Kupka et al. 1999). For a few elements 
it was necessary to make corrections to the atomic line 
data, and these corrections were made during the final 
steps of abundance analysis (see sections below for indi- 
vidual elements for details). The spectrum synthesis code 
SYNTH (Piskunov 1992) was used in all synthetic spec- 
trum calculations. Using these calculations we estimated 
rotational velocities for both stars, which were found to 
be wsini=9±2 kms" 1 for HD 32115 and usini=17±2 
kms -1 for HD 37594. While for the most part rotational 
plus instrumental broadening provides adequate fits to un- 
blended spectral lines in HD 32115, this is not the case for 
HD 37594. Single spectral lines in this star often have V- 
shaped profiles rather than the usual U-shaped profiles 
produced by the rotational broadening. Similar line pro- 
files are observed in pulsating stars (eg. the study of the 
6 Set- type star FG Vir by Mittermayer 2001). This effect 
is marginal, and needs a further study with higher spectral 
resolution. 

The final abundance analysis was done by three dif- 
ferent methods: classical analysis using equivalent width 
data, spectral synthesis for strong and partially blended 
lines, and a semi-automatic procedure which uses both 
equivalent widths and spectrum synthesis techniques. The 
results of the first two methods were combined. We found 
no systematic difference between equivalent widths mea- 
sured in the Solar Flux Atlas (Kurucz et al. 1984) and in 
a day-sky spectra obtained at the same dates and with 
the same spectrograph that was used for the taking the 
spectra of the program stars. The range of the errors in 
equivalent width measurements is 2 5 mA and depends 
on S/N ratio (position of the line on the echelle order), and 
blending effects. For most of the unblended lines this error 
does not exceed 2-3 mA. We also compared two indepen- 



dent sets of equivalent width measurements in HD 32115 
spectrum made by I. Bikmaev with the Gaussian approx- 
imation and by H. Bruntt with the Voigt approximation. 
For 122 lines the mean difference between two sets is 1.2 
mA with a standard deviation 2.8 mA. 



5.1. Classical abundance analysis 



Comparing the synthetic and the observed spectra we se- 
lected spectral lines which are free of substantial blending 
effects and measured their equivalent widths. The final 
abundances were calculated with Kurucz's WIDTH9 code 
modified by V. Tsymbal for VALD-format input files. For 
a few elements, like Ban which has mostly very strong 
lines, or Euli and Ndni which have only 1-2 partially 
blended lines in the whole spectral region, abundances 
were derived with synthetic spectrum calculations. Final 
abundances for 31 chemical elements (41 ions) are col- 
lected in Table 3 (columns 2 and 7 for HD 32115 and 
HD 37594 respectively). The standard deviation is given 
in units of the last digit in parentheses. Solar abundances 
in \og(N/N to t) = \og(N/N H ) - 0.04 (correction due to 
He) are taken from Grevesse & Sauval (1998) except C, 
N, O, Mg, Si, and Ba. For these elements, except Ba, 
the most recent data by Holweger (2001) are given. We 
adopted the solar Ba abundance derived by Mashonkina 
& Gehren (2000). NLTE corrections were applied for the 
species marked by asterisks (see details below) . Columns 6 
and 11 contain the stellar abundance relative to the sun. 
We also calculated line-by-line differences between stel- 
lar and solar abundances, which show a close agreement 
with the values from columns 6 and 11. All individual 
line abundances are presented in Table 4 which is avail- 
able in electronic form only. This table may be retrieved 
from http://www.inasan.rssi.ru/~ryabchik. NLTE correc- 
tions are also included in Table 4. These are probably the 
most complete abundance data for single A-type stars to 
date. 
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5.2. Semi-automatic abundance analysis 

We have written software called vwa for semi-automatic 
abundance analysis which is described in Bruntt ct al. 
(2002). The code runs though the list of atomic line data 
extracted from vald (Kupka et al. 1999) for the target 
star, and selects lines that suffer from the least amount of 
blending. For a region around each selected line the syn- 
thetic spectrum is calculated for different abundances of 
the element forming the line until the equivalent widths of 
the observed and synthetic spectrum match. The compu- 
tation of the synthetic spectrum is done with SYNTH ver- 
sion 2.5 (Valenti & Piskunov 1996). The code needs the 
line data from vald and the adopted atmosphere model 
as input. When vwa has fitted the selected lines they are 
inspected visually - and the fit is either accepted, rejected 
or must be improved manually. 

While classical abundance analysis can be performed 
relatively fast for a star like HD 32115 which has low 
vsmi the advantage of using vwa for abundance anal- 
ysis is demonstated for HD 37594. Since this star has a 
moderate value of vsini one needs to select lines with 
some amount of blending to be able to derive abundances 
of as many elements as possible. We first used VWA where 
we only selected unblended Fe lines. We then adjusted the 
microturbulence and recalculated abundances with VWA 
until the Fe I abundance and equivalent width of the lines 
did not correlate. We thus found £ t = 2.4 ± 0.3 kms -1 . 
Since the abundance of Fe was now considered to be ac- 
curate, we allowed vwa to select lines of other elements, 
which were mild blends with Fe. By this method we were 
able to use lines that could not be used in the classical 
abundance analysis which is strictly limited to unblended 
lines. 

When comparing the results of the abundance analy- 
sis of the two methods in Table 3 we see that far fewer 
lines are chosen by vwa but the derived abundances all 
agree within the error estimates. Hence, vwa is usable for 
carrying out fast and reliable abundance analysis. 

5.3. Errors in abundance determinations 

Our derived abundances contain errors from a number 
of sources of uncertainty, including the adopted atomic 
data, most importantly the oscillator strengths, the equiv- 
alent width measurements, the adopted atmospheric pa- 
rameters and model, and the observed spectra. For some 
sources we can directly estimate the associated error. The 
estimated uncertainty in equivalent width measurements 
is ±2 mA, which results in 0.20-0.09 dex error in the abun- 
dance for lines with equivalent widths 5 to 10 mA, 0.05 
dex for W x =20 mA, and less than 0.03 dex for W\>40 
mA. An error of 100 K in the effective temperature deter- 
mination results in a systematic error of 0.05 dex on the 
the abundance of the most temperature-sensitive atomic 
species, e.g. Ti I, V I. For cases where a suitably large num- 
ber of lines are used, the standard deviations quoted in 
Table 3 estimate the total random error associated with 



all sources of uncertainty. However, for species where only 
a small number of lines are available this cannot be ex- 
pected to give any indication of this error. One can see 
from Table 3 that the standard deviations for most ele- 
ments with a sufficient number of spectral lines do not 
exceed 0.15 dex. Considering this as indicative for all el- 
ements, and taking into account also possible systematic 
errors we estimate the total error at about 0.2 dex for the 
species with 1-3 spectral lines. 

Formally, a gravity decrease by 0.2 dex may provide 
ionization balance for all iron peak elements, since abun- 
dances obtained from the ionized lines will be decreased 
by 0.06-0.09 dex, while abundances from the neutral lines 
would not be changed. CNO abundances would be reduced 
by 0.05 dex. However, we stress that the effective gravity 
for both stars is accurately defined by parallax measure- 
ments. 

Since most solar photosphcric abundances are derived 
with the semi-empirical model of Holweger & Muller 
(1974), it is important to check that our use of the- 
oretical Kurucz models does not itself introduce any 
significant differences relative to the solar abundances. 
To test this we compared solar Fe abundances derived 
from Fen lines in the solar spectrum employing the 
Holweger-Mullcr solar atmospheric model (Holweger & 
Muller 1974) with those found using the solar model com- 
puted with ATLAS9 and the same convection treatment as 
for our stars (see Gardiner et al. 1999). The abundances 
were calculated from the best 13 Fen lines (Schnabel, 
Kock & Holweger 1999) with oscillator strengths from 
their paper and from the vald database. The results 
are: logA^ Fc =7.42±0.09 (SKH log 5 / values, HM model, 
weighted mean)), log N Fc =7.46±0.11 (SKH \oggf values, 
our model, unweighted mean), log N-p c =7.48±0.10 (vald 
log gf values, our model, unweighted mean). With our so- 
lar model the line-by-line difference (HD 32115— Sun) for 
31 common Fen lines is -0.02±0.10 and agrees within the 
error limits with the [Fe/H] value presented in Table 3. 

5.4. CNO elements 

Oscillator strengths for C and N were taken from Hibbert 
et al. (1993) and Hibbert et al. (1991) respectively. The 
same data were used in reevaluation of the solar abun- 
dances (Holweger 2001). For oxygen, we employed the os- 
cillator strengths supplied by vald which originate from 
the NIST compilation (Wiese et al. 1996). We did not use 
the IR O I triplet because of the possibly strong NLTE 
effects. Expected NLTE corrections are -0.1 dex for C 
(see Sturenburg, 1993, Rentzsch-Holm 1996, Paunzcn ct 
al. 1999) and -0.01 - -0.03 dex for O (Takcda 1997). Both 
stars show similar small deficiencies of C and N. 

5.5. A/a to K 

These elements are known to be subject to NLTE effects. 
We also note that the oscillator strengths for Sil lines in 
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Table 3. Abundances in HD 32115 and HD 37594 based on n measured lines with the standard deviation in units of the last 
two digits in parentheses. NLTE corrections are applied for the species marked by asterisks. The letters 'H' and 'MG' in the 
last column mean that solar abundances are taken from Holweger (2001) and Mashonkina & Gehren (2000); otherwise standard 
solar composition as reported by Grevesse & Sauval (1998) is used. 



Ion 


HD 32115 


HD 37954 


Sun 




Classical an 


alysis 


Semi-automatic 




Classical a 


lalysis 


Semi-automatic 








log(JV/JV to t) 


n 


\og{N/N t ot) 


n 


[N/Nh] 


log(A7iV t0 i 


n 


log{N/N to t) 


n 


[N/Nh] 


\og(N/N t ot) 


C I 


-3.62(17) 


35 


-3.52(24) 


5 


-0.17 


-3.73(15) 


27 


-3.74(05) 


4 


-0.28 


-3.45H 


N I 


-4.19(08) 


11 


-4.22(04) 


2 


-0.08 


-4.27(12) 


4 






-0.16 


-4.11H 


O I 


-3.23(10) 


12 


-3.25(04) 


2 


0.07 


-3.36(12) 


5 


-3.27(12) 


5 


-0.06 


-3.30H 


Na I* 


-5.86(12) 


8 






-0.15 


-6.05(21) 


8 






-0.34 


-5.71 


Mgl 


-4.59(12) 


9 


-4.33(05) 


4 


-0.09 


-4.89(10) 


8 


-4.68(03) 


4 


-0.39 


-4.50H 


Mg II 


-4.43(11) 


4 


-4.46 


1 


0.07 


-4.79(19) 


3 


-4.65 


1 


-0.29 


-4.50H 


Al I 


-5.69(15) 


4 






-0.12 


-6.15(06) 


2 






-0.58 


-5.57 


Si I 


-4.53(13) 


39 


-4.49(12) 


11 


-0.03 


-4.75(08) 


26 


-4.76(13) 


11 


-0.25 


-4.50H 


Si II 


-4.47(12) 


6 


-4.53 


1 


0.03 


-4.59(04) 


3 


-4.72(10) 


2 


-0.09 


-4.50H 


P I 


-6.42(13) 


2 






0.17 












-6.59 


S I 


-4.81(15) 


19 


-4.92(09) 


5 


-0.10 


-5.07(10) 


8 


-5.06(01) 


2 


-0.36 


-4.71 


CI I 


-6.76 


1 






-0.22 












-6.54 


K I* 


-7.05 


1 






-0.13 


-7.32 


1 






-0.40 


-6.92 


Cal 


-5.63(15) 


21 


-5.60(13) 


13 


0.05 


-5.89(22) 


26 


-5.99(27) 


16 


-0.21 


-5.68 


Ca II* 


-5.69(11) 


7 


-5.58 


1 


-0.01 


-5.89(04) 


4 


-5.88 


1 


-0.21 


-5.68 


Sc II 


-8.89(10) 


11 


-8.96(20) 


4 


-0.02 


-9.04(15) 


6 


-9.13 


3 


-0.17 


-8.87 


Ti I 


-7.20(14) 


21 


-7.18(04) 


4 


-0.18 


-7.43(14) 


3 


-7.42(03) 


2 


-0.41 


-7.02 


Ti II 


-7.03(15) 


27 


-6.92(20) 


13 


-0.01 


-7.23(18) 


19 


-7.38(12) 


8 


-0.21 


-7.02 


V I 


-8.24(25) 


3 






-0.20 












-8.04 


V II 


-8.17(08) 


5 


-7.95 


1 


-0.13 


-8.13(25) 


2 






-0.09 


-8.04 


Cr I 


-6.40(12) 


39 


-6.40(13) 


7 


-0.03 


-6.71(10) 


11 


-6.75(06) 


5 


-0.34 


-6.37 


Cr II 


-6.27(09) 


20 


-6.22(12) 


7 


0.10 


-6.63(09) 


17 


-6.55(06) 


6 


-0.26 


-6.37 


Mn I 


-6.87(17) 


20 


-6.76(37) 


6 


-0.22 


-7.14(16) 


9 


-7.04(13) 


4 


-0.49 


-6.65 


Mn II 


-6.72(11) 


2 






-0.07 












-6.65 


Fe I 


-4.60(12) 


235 


-4.59(11) 


109 


-0.06 


-4.88(13) 


183 


-4.89(12) 


38 


-0.34 


-4.54 


Fe II 


-4.48(14) 


60 


-4.54(11) 


23 


0.06 


-4.87(11) 


25 


-4.76(11) 


4 


-0.33 


-4.54 


Co I 


-7.16(10) 


5 


-7.16(00) 


2 


-0.04 


-7.43(03) 


3 






-0.31 


-7.12 


Ni I 


-5.91(10) 


59 


-5.94(13) 


17 


-0.12 


-6.24(16) 


27 


-6.18(03) 


4 


-0.45 


-5.79 


Ni II 


-5.82 


1 






-0.03 












-5.79 


Cu I 


-7.82(26) 


5 






0.01 


-8.32(06) 


2 






-0.49 


-7.83 


Zn I 


-7.77(16) 


3 






-0.33 


-8.07(05) 


3 






-0.63 


-7.44 


Sr II* 


-9.09(22) 


3 






-0.02 


-9.19(07) 


3 






-0.12 


-9.07 


Y II 


-9.68(15) 


8 


-9.74(08) 


3 


0.12 


-9.94(08) 


8 


-9.95(07) 


3 


-0.14 


-9.80 


Zr II 


-9.38(16) 


7 






0.06 


-9.42(20) 


5 






0.02 


-9.44 


Ba II* 


-9.64(12) 


6 


-9.57(05) 


3 


0.19 


-9.87(28) 


4 


-9.70(19) 


4 


-0.04 


-9.83MG 


La II 


-10.82(08) 


4 






0.05 


-11.15(13) 


3 






-0.28 


-10.87 


Ce II 


-10.41(19) 


7 






0.05 


-10.57(16) 


6 






-0.11 


-10.46 


Nd II 


-10.49(25) 


11 






0.05 


-10.63(01) 


2 






-0.09 


-10.54 


Nd III 


-10.40(07) 


2 






0.14 












-10.54 


Sm II 


-11.06 


1 






-0.03 


-11.38 


1 






-0.35 


-11.03 


Eu II 


-11.63 


1 






-0.10 


-11.80 


1 






-0.27 


-11.53 



vald are not accurate, many of them coming from Kurucz 
& Peytremann (1975). Therefore we recalculated oscilla- 
tor strengths for most of the Si I lines by fitting them to 
the Solar Flux Atlas (Kurucz et al. 1984) adopting the 
solar abundance in Table 3. NLTE corrections were calcu- 
lated for Nai, and K I following atomic models and proce- 
dures developed by Mashonkina et al. (2000) for Na, and 
Ivanova & Shimansky (2000) for K. Oscillator strength 
data for All were taken from Baumuller & Gehren (1996). 



According to this paper NLTE corrections for the All 
lines used in our abundance analysis are negligible. For 
Mgl lines AA 4703, 4730, 5528, 5711 oscillator strengths 
were taken from Jonsson et al. (1984) and from Froese 
Fischer (1975). In both stars the Mg abundance was de- 
rived by the spectral synthesis technique with the Stark 
and Van der Waals damping constants giving the best fit 
for the same Mg I lines in the spectra of Procyon and the 
Sun (Fuhrmann et al. 1997, Zhao et al. 1998). Individual 
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LTE and NLTE abundances are given in Tabic 4. The 
largest NLTE corrections were found for the Kl reso- 
nance line A 7698.97. They are -0.61 (HD 32115) and -0.41 
(HD 37594). Unfortunately, another Kl resonance line at 
A 7664.91 falls in one of the gaps in our cchcllc spectra. 
We could measure the strongest Cli line A 8375.95 in the 
spectrum of HD 32115. The solar photospheric CI abun- 
dance is very uncertain, log(C7/7V tot )=-6.54±0.3. The CI 
abundance in HD 32115 was obtained from just one weak 
line but agrees perfectly with the CI abundance found in 
meteorites (Grevesse & Sauval 1998). Two strongest in- 
frared Pi lines AA 9593.50, 9796.83 measured in the latest 
spectrum of HD 32115 were used to derive phosphorus 
abundance. 

5.6. Ca to Ni 

The iron-group elements usually represent the metallicity 
of a star. With a few exceptions we were able to obtain 
abundances from the lines of neutral and singly-ionized 
species for all elements of this group, thus checking the 
ionization balance in the atmospheres. The abundances 
obtained by the neutral and ionized Fe-group elements 
differ by «0.13 dex on average. In both our stars the Fe- 
group elements in the ionized state are dominating and 
therefore their lines are less sensitive to temperature inac- 
curacy and NLTE effects and hence should provide better 
abundance estimates. Among iron-group elements NLTE 
calculations were carried out for Ca (see below) and Fe 
(Rentzsch-Holm 1996). An expected NLTE correction to 
the abundance obtained from Fei lines is «+0.1 dex for 
T e ff«7300 K. Applying this correction as a first approxi- 
mation to the abundances obtained from the lines of neu- 
tral elements a reasonable ionization balance is obtained. 

Five Fe n lines correspondning to 4d-4f transitions with 
an upper energy level at around 10.5 eV were measured 
in the HD 32115 spectrum. For these lines only theoreti- 
cal oscillator strengths are available. Both the log gf val- 
ues from the Kurucz line list (included in vald) (Kurucz 
1993) and the more accurate sets of calculated \oggf- 
values from Raassen & Uylings (1998) agree to within 0.07 
dex. Abundances obtained from these lines do not deviate 
by more than 0.1-0.2 dex from the mean Fe abundance 
and thus supports the oscillator strengths for these high 
excitation Fe II lines to this level of accuracy. 

Fe abundance analysis in HD 32115 is based on about 
300 lines well spread through the whole spectral region of 
4100-9200 A. We did not find any clear dependence of the 
individual abundances on the wavelength (Fig. 2) which 
might be expected in the primary star if the contribution 
from the secondary is not negligible. 

For the two strongest Call lines AA 8912, 8927 we in- 
creased logg/ by 0.3 dex to fit line profiles in the so- 
lar spectrum. Stark damping constants for these lines, 
which are important, are taken from Dimitrijevic & Sahal- 
Brechot (1993). We also carried out NLTE calculations for 
Ca II lines using a model of the atom developed by Ivanova 



et al. (2002). NLTE corrections are always negative and 
of the order of -0.15 dex for AA 8248, 8912, 8927 lines in 
both stars, while they are negligible for weak Can lines. 



-3.5 



-4.0 - 




-5.5 r i 

4000 5000 6000 7000 8000 9000 10000 
Wavelength 

Fig. 2. The dependence of the individual abundances obtained 
from Fei lines in HD 32115 on wavelength. 



5. 7. Cu and Zn 

Abundances of these elements were obtained by spectral 
synthesis taking into account hyperfine splitting (hfs). Hfs 
constants were extracted from Biehl (1976). Solar syn- 
thetic spectrum calculations showed that for Zn I lines we 
may neglect the hfs effect, while it is strong for Cu I lines in 
the solar spectrum. In both our stars Cul lines are weak, 
and the hfs effect for Cul and Znl lines is not significant. 

Zn shows the largest deviation from the solar abun- 
dance in HD 32115, and the largest deficiency in 
HD 37594. To our knowledge, there are no NLTE calcula- 
tions for Cul or Znl. With ionization energies \ = 7.72 eV 
for Cul and 9.39 eV for Znl these atoms are minor species 
similar to Mgi, All, Fei in the stellar atmospheres typical 
for investigated stars. Based on our experience in NLTE 
analyses and taking into account NLTE calculations for 
atoms with similar atomic parameters and term structures 
(Sturenburg 1993, Rentzsch-Holm 1996) Cul and Znl are 
expected to be over-ionized. In this case low atomic levels 
are underpopulated compared with LTE, and the lower 
the excitation energy is the greater the departures from 
LTE are. Most probably, NLTE effects will weaken the 
investigated Cu I and Zn I lines and NLTE abundance cor- 
rections are expected to be positive. Usually NLTE correc- 
tions are increasing with the line intensity. The observed 
equivalent widths of Znl lines are 2-3 times larger than 
for Cul lines, therefore we may also expect larger positive 
NLTE corrections for Zn. 
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5.8. Sr, Y, Zr, Ba 

NLTE corrections for Sr n lines were calculated using the 
model atom developed by Belyakova et al. (1997, 1999). 
They do not exceed 0.1 dex. Hfs effects may be impor- 
tant for Yn lines, but there are no data available. Zr 
abundances are rather uncertain, because reliable oscilla- 
tor strengths are only known for a small number of spec- 
tral lines. We could not use the solar spectrum fitting for 
Zrii lines measured in our stars due to strong molecular 
contributions in the corresponding spectral regions. 

NLTE corrections to the Ba abundance were calculated 
according to Mashonkina et al. (1999). They are large, up 
to 0.18 dex, for the two strongest Ban lines AA 6141, 6496. 
It should be mentioned that the solar Ba abundance used 
by us was derived by Mashonkina & Gehren (2000) from 
the same lines as in our analysis, and it is closer to abun- 
dance of Ba found in meteorites than any previously pub- 
lished solar Ba abundance values (see Grevesse & Sauval 
1998). Stark damping constants for Ban lines were taken 
from Dimitrijcvic & Sahal-Brcchot (1996), while for reso- 
nance Sr II lines the corresponding values were taken from 
the NIST compilation (Konjevic et al. (1984)). 

5.9. La to Eu 

New experimental oscillator strengths and hfs constants 
were used for Laii (Lawler et al. 2001a), Ceil (Palmcri 
et al. 2000, Zhang et al. 2001), and Euli lines (Lawler 
et al. 2001b). Two weak lines of Ndm AA 5102, 5295, 
which are the strongest lines in the Nd in spectrum, were 
identified in HD 32115. They are partially blended so 
abundance estimates were made by spectral synthesis. 
Oscillator strengths and other atomic parameters in VALD 
are taken from Cowley & Bord (1998) and from Bord 
(2000). Both lines provide Nd abundances close to that ob- 
tained from the analysis of Nd II lines, thus giving strong 
support for the Nd in oscillator strength calculations. 
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Fig. 3. The observed relative abundances in HD 32115 (filled 
circles) and in HD 37594 (open circles). 
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Fig. 4. A comparison between Ha line profiles observed in 
HD 32115 (thick full line) and in roAp star HD24712 (dashed 
line), and calculated for the adopted 7250-ODF model (thin 
full line). 



6. Discussion 

The abundances found for HD 32115 and HD 37594 rela- 
tive to the sun are shown in Fig. 3 by filled and open circles 
respectively. Within the typical error limits of ±0.15 dex 
HD 32115 is a solar abundance star (a mean metallicity 
[M] = -0.04 ± 0.11). Therefore this star may be used 
as a chemical standard in studies of cool peculiar stars. 
We also get a standard for further investigation of the hy- 
drogen wing-to-core anomaly found in cool Ap stars by 
Cowley et al. (2001). Fig. 4 shows a comparison between 
Ha line profiles in HD 32115 and in one of the pulsat- 
ing Ap (roAp) stars, HD 24712, with the same effective 
temperature. It is evident from this comparison that the 
Ap-star anomaly occurs in that part of the hydrogen line 
core which is not reproduced by our current models. 

HD 37594 is slightly metal-deficient. Its mean metal- 
licity [M] = -0.26 ± 0.16 is close to the value [M] = 
— 0.15±0.04 obtained from photometric calibrations. This 



error in the metallicity obtained from photometric cali- 
brations is based only on the errors in the photometric 
indices. No systematic errors of the calibration itself are 
included, which may be up to 0.1. Further support for 
the high accuracy of the metallicities derived from the 
Stromgren photometry for normal stars is provided by 
the results of abundance analysis of 28 And = HD 2628 
(Adelman et al. 2000) and o Boo A=HD 128167 (Adelman 
et al. 1997). [M]=-0.16±0.29 and -0.29±0.37 were derived 
for these stars from abundances versus -0.13 and -0.35 
from photometric calibration. The abundance pattern in 
HD 37594 is similar to that in 28 And (T off =7250 K, 
log g=4.2, £ t =2.3 kms -1 , vs'mi=9 kms -1 - see Adelman 
et al. 2000) which belongs to the S Set group. The smaller 
scatter derived in the present paper is explained by more 
accurate abundance determinations with NLTE and hfs 
effects taken into account as well as more careful choice of 
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the unblended lines which is possible in the red spectral 
region. 

Abundances obtained for HD 32115 and HD 37594 do 
not follow the predicted abundance pattern for Z <28 
from consistent stellar evolution models calculated with 
radiative forces, opacities and diffusion (Turcotte et al. 
1998). Models for 1.5 M/Mq predict that elements with 
5 < Z < 20 will be underabundant by 0.5 dex relative to 
Fe, and that iron-peak elements will be generally slightly 
overabundant even at early evolutionary phases on the 
MS. This was not observed in either of the stars analysed 
here. 

To confirm the derived abundances we have used the 
software vwa (Bruntt et al. 2002) which was developed 
to be able to make fast semi-automatic abundance anal- 
ysis. We find the same abundances as the more careful 
classical approach. At present vwa is being applied to the 
study of the primary target candidate stars for the aster- 
oseismology missions COROT and R0mer (Bruntt et al. 
2002). Several of these stars have moderate or high vsini 
and the amount of work needed to carry out the analysis 
using vwa is reduced substantially. 
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